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We examine the correlation between supernova host galaxy properties and 
their residuals on the Hubble diagram. We use supernovae discovered during 
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the Sloan Digital Sky Survey II - Supernova Survey, and focus on objects at 
a redshift oi z < 0.15, where the selection effects of the survey are known to 
yield a complete Type la supernova sample. To minimize the bias in our anal- 
ysis with respect to measured host-galaxy properties, spectra were obtained for 
nearly all hosts, spanning a range in magnitude of —23 < Mr < —17. In contrast 
to previous works that use photometric estimates of host mass as a proxy for 
global metallicity, we analyze host-galaxy spectra to obtain gas-phase metallici- 
ties and star-formation rates from host galaxies with active star formation. From 
a final sample of ~ 40 emission-line galaxies, we find that light-curve corrected 
Type la supernovae are ~ 0.1 magnitudes brighter in high- metallicity hosts than 
in low- metallicity hosts. We also find a significant (> 3cr) correlation between 
the Hubble residuals of Type la supernovae and the specific star-formation rate 
of the host galaxy. We comment on the importance of supernova/host-galaxy 
correlations as a source of systematic bias in future deep supernova surveys. 

Subject headings: cosmology: observations — supernovae: general — surveys 



Introduction 



The utility of Type la supernovae (SNe la) to the study of cosmology rests upon the 
ability to calibrate the difference in absolute magnitude between events, allowing for an 
accurate determination of the distance-redshift relationship. Variations in the intrinsic lu- 
minosity of ea c h event can be measured by the widt h and color of the observed light curve 
( jPhillipsI Il993t iHamuy et al.l Il996l : iRiess et al.lll996l). A fter applying these relations using 
light-curve fitting codes such as salt2 (IGuv et al.lboOTI ) or mlcs2k2 (jjha et al.lboOTl ). the 
resulting derived distance moduli display an intrinsic scatter about the best-fit cosmology 
equivalent to ~ 7% in distance. The difference between the measured distance modulus of a 
SN and the best-fit cosmology is known as the Hubble Residual (HR). 

The absolute magnitude of SNe la have additionally been shown to depend on their 
environment; SN e la in late-type g a laxies are intririsically more luminous than those in 
early-type hosts ( iHamuy et al.l Il996l : [Gallagher et al.ll2005l). On e possible explanation for 



this difference is the progenitor metallicity. iTimmes et al.l (120031 ) showed, both analytically 
and through modelling, that a metal-rich progenitor produces less ^^Ni, and as such is less 
luminous, than a metal-poor SN la progenitor. If this difference in luminosity as a function 
of metallicity is accompanied by the same changes in color and light-curve stretch that 
normally occur in SNe la , then this environ mental effect would not produce biased distance 
measurements. However, iKasen et al.l (120091 ) have recently shown that light-curve corrections 
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tend to overcompensate for the metallicity effect. Specifically, they show that for a given 
quantity of synthesized ^^Ni, multi-dimensional models predict that a higher metallicity 
progenitor will produce a narrower light curve. This effect means that after standard light- 
curve corrections are applied, a high-met allicity progenitor would appear to be overluminous. 

Several groups have recently shown that the host galaxy environment does in fact 
correlate with the HR, in the sense that more massive galaxies host overluminous (for 



their light-curve shap e) SNe la flGallagher et al.ll2008l : iKellv et al.ll20ld : ISullivan et al.ll2010 



Lampeitl et al.ll2010l ). Under the assurnption that these galaxies all follow the mass-metallicity 
relationship deri ved inlXremoiiti et al" ( 2004 , hereafter T04), where high-mass galaxies have 
high metallicity, ISullivan et al.l (I20101) showed th at their results are consistent with the metal- 
licity /luminosity relation in iKasen et al.l ( 120091 ). 



Knowing the true cause of the HR correlation with host galaxy properties is key for fu- 
ture deep SN surv eys such as the Dark Energy Survcy EI (DES) and the Large Synoptic Survey 
Telescope (LSST; ILSST Science Collaboration! 120091 ). While it is clear that correcting for 
this relationship would lead to decreased scatter in HRs, care must be taken to consider 
the change in galaxy properties with redshift. If progenitor metallicity truly is the cause of 
the observed correlation, and host mass is used as a proxy for t he metallicity, then e volu- 
tion in the mass-metallicity relationship as a function of redshift (IMaiolino et al.ll2008l ) will 
bias the Hubble Diagram and our measurement of cosmological parameters. The correlation 
might instead be due to the relative prevelance of progenitor channels and thus depend on 



star-fo rmation rate (SFR), which will have a different redshift dependence. iGallagher et al. 



( I2OO8I ) showed a correlation exists between stellar metallicities from passive galaxies and SN 
la HRs at low redshifts {z < 0.05), but no correlation between measu red emission-line meta l- 
licities from star- forming galaxies (SFG) has yet been observed (see [Gallagher et al.ll2005l ). 
Furthermo re, other factors m ay contribute part of the intrinsic scatter, such as explosion 
geometry (IMaeda et al.l 1201 II ). without being directly tied to host-galaxy properties. 



In what follows we determine the relationship between HRs and both host galaxy metal- 
licity and SFRs in a nearly unbiased sample of spectroscopically-observed host galaxies of 
SNe la discovered as part of the SDSS-II Supernova Survey. Section [2] describes our ob- 
servations and the construction of our sample. We discuss our method for measuring the 
gas-phase metallicity and specific SFR (SFR per unit stellar mass; sSFR) of star-forming 
host galaxies of SNe la in Section [31 Results are presented in Section HI and the implications 
of our findings for SN la cosmology are detailed in Section [51 



^http: / /www. darkenergysurvey.org/ 
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Observations 



We apply our analysis to host galaxies of the SNe la disc overed as part of the Sloan 
Digital Sky Survey-II Supernova Survey (hereafter SDSS-SNS: iFrieman et al.ll2008[ ). The 
SDSS-SNS repeate dly surveyed the 300 square degree Southern Equatorial Stripe (desig- 
nated stripe 82; see IStoughton et al.ll2002l ) during the Fall seasons (September 1 - November 
30) of 2005-2007 usi ng the dedicated 2.5 meter SDSS telescope at Apache Point Obser- 
vatory, New Mexico (IGunn et al.l |2006| ) . Each photometric obse rvation consists of ne arly 
simultaneous 55 second exposures in ea ch of the five uqri z filters ( jFukugita et al.lll996l ) us- 
ing the w ide-field SDSS CCD camera (IGunn et al.lll998l ). High quality light curves were 
obtained (IHoltzman et al.ll2008l) on a photometric system calib rated to an uncert ainty of 1% 
(jivezic et al.l 12004 ) . For a technical summary of the SDSS, see lYork et al.l ( 120001 ). 



The SDSS-SNS spectroscopically confirmed 504 SNe la over the duration of its three 
year survey. Details of the spectros copic observing c ampaign and the criteria by which 
targets were selected can be found in ISako et al.l (120081 ). An additional 210 transients with 
identifiable hosts have been designated photometrically probable SNe la. These objects 
did not have a spectrum taken during the survey, usually because they either (a) were not 
well separated from their host, or (b) were below the magnitude limit for our foUowup 
spectroscopy. However, these objects are highly likely to be SNe la as opposed to any of the 
core-collapse types based on their multicolor light curves and the redshifts obtained from 
spectra of their host galaxies. The u pdated photometr ic classifier code, containing more core- 
collapse templates and described in ISako et al.l (1201 if ), was designed for the S DSS-SNS and 



has b een shown to be accurate at determining SNe la with low contamination ( iKessler et al. 



20101 ). Where necessary we will refer to the spectroscopically-confirmed sample as Spec la 
SNe and the photometrically-probable sample with identified hosts as Phot la SNe. 



We derive distance moduli from our sample of SNe la with the SNANA code (IKessler et al 



2009bl ). using both the mlcs2k2 and salt2 light-curve fitters. Light-curve quality cuts 
applied are the sarn e as those in Section 4 of the SDSS-SNS first-year cosmology paper 
( IKessler et al.ll2009al ). with the one exception being a more stringent requirement of at least 
one measurement at two days or more before peak brightness in the rest-frame 5-band ac- 
cording to mlcs2k2 (Trcst < —2). These cuts remove SNe la that have low signal-to-noise 
measurements, insufficient temporal coverage, and peculiar light-curve sh apes. Distanc e 



moduli from SALT2 light-curve fits are computed using the code 'SALT2mu' (IMarriner et al. 



20111). The values used for the a and /3 paramet ers, which are the corrections for stretch 
and color, are determined in iMarriner et al.l ( 120111 ) as those best-fit by the SDSS-SNS data, 
independent from cosmology. For mlcs2k2, the reddening law Ry = 2.03 is used, as is the 
default SNANA prior on Ay of exp (— ^4^/0.3). Hubble Residuals are determined by sub- 
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tracting the distance modulus of the assumed cosmology at the redshift of the host galaxy 
from that of the SN (HR = /xrn— /x^). W e use the best-fit ACDM cosmology to the SDSS-only 



SN sample from iKessler et al.l (l2009al ). D.m = 0.274 and = 0.735. We assume a Hubble 
Constant of Hq = 70 km sec~^ Mpc~^ throughout this paper, though we note that the choice 
of Hubble Constant is irrelevant to our results, as Hq is degenerate with the fiducial SN la 
absolute luminosity. 

A point of emphasis in our analysis is to ensure that selection biases in our data set are 
minimized to the greatest extent possible. To this end we have limited our study to SNe la 
at redshifts z < 0.15, where it has been demonstrated that the selectio n efficiencies of the 



surve y are such that the Spec la plus Phot la SNe sample is complete ( iDilday et al.l 12008 



20101 ). Thus our results will not be biased by containing an over-abundance of intrinsically 
over-luminous SNe la. The galaxy spectra analyzed in this paper are hosts of the 140 Spec 
la and 7 Phot la SNe at 2; < 0.15, of which 77 Spec la and 3 Phot la SNe pass our light-curve 
quality cuts. 

A concerted effort is being undertaken to obtain galaxy-only spectra for all Spec la and 
Phot la SNe discovered in the SDSS-SNS. For the low-redshift sample analyzed in this paper, 
this campaign is essentially complete; we have host-galaxy spectra for 93% percent of the low- 
z sample. Over half of these SNe la hav e host-galaxy spectra f rom the SDSS Legacy Survey 



(hereafter referred to as 'SDSS spectra'; lAbazajian et al.ll2009[ ). and were obtained from the 



SDSS Data Archive Server ( DAS). In collaboration with the SDSS-III Baryon Oscillation 



Spectroscopic Survey (BOSS; lEisenstein et al.ll201ll ) most of the remaining unobserved SNe 



la hosts with rrir < 20.5 were targeted, including over a third of the low-2; sample. The 
faintest host galaxies were observed with the 8m Gemini-South Telescope and are vital for 
the completeness of our sample, as the least luminous galaxies tend to have the lowest 
metallicities. The remainder of the host-galaxy spectra used in this analysis were obtained 
with the 3.5m Astrophysical Research Consortium Telescope at Apache Point Observatory 
(APO) and the 3.6m New Technologies Telescope (NTT), which observed some hosts lacking 
SDSS spectra before the BOSS and GEMINI observations were begun. 

The fact that we begin with a nearly complete sample of host-galaxy spectra is crucial, 
as this means we are not biased toward more luminous hosts. The parameters that we are 
attempting to directly measure (metallicity, SFR) are correlated with the absolute magnitude 
of a galaxy, as demonstrated in T04: metal-poor galaxies tend to be faint, metal-rich galaxies 
bright. It is thus necessary to analyze hosts with a wide range of absolute magnitudes if we 
intend to measure a wide range of progenitor metallicities. We demonstrate in Section 14.31 
that the final sample we use in this paper, after all of our data cuts are apphed, still satisfies 
this criterion, though we are left with a much reduced data set. 
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2.1. Spectroscopic Data Processing 

The primary focus of this section is to describe the reduction of data obtained at the 
G emini Observat o ry. In formation pertaining to the reduction of N TT spectra can be fo und 



m 



Ostman et al.l f l2010l ). while the APO procedure is described in IZheng et al.l (|2008[ ). A 
forthcoming paper (Olmstead et al., in prep) will describe the data pipeline from BOSS. All 
spectra from the SPSS spectroscopic survey used in this paper are from data release 7 (DR7; 



Abazajian et al.l 120091 ). 



The faintest identified host galaxies of SNe la at z < 0.15 (determined from SN spectra) 
were observed with the Gemini Multi-Object Spectrograph (GMOS) at the Gemini-South 
observatory. We were awarded 23 hours of time during semester 2008B (program GS-2008B- 
Q-38), during which 17 host spectra were obtained. We took spectra in longslit mode with 
a 1.0 arcsec slit, no filter, and used the B600 grating. The approximate observed wavelength 
coverage was 3800-6700 A. Three exposures of duration 500-800 seconds each were taken for 
central wavelengths of both 520 and 525 nm, where the offset mitigates the effect of bad 
detector pixels in our reduced data. The data were binned 2 times in the spatial direction 
and 2 times in the spectral direction for a wavelength dispersion of 0.9 A pixel"^. The B600 
grating has a spectral resolution of 1688 at its blaze wavelength of ~ 4610A when used with 
a 0.5 arcsec slit. 

We performed most of the data reduction using version 1.10 of the Gemini IRAfI 
package. The raw two-dimensional (2D) spectra are in multi-extension FITS (MEF) format; 
we propagate the variance (VAR) and data quality (DQ) planes throughout our reduction. 
We created mean bias frames by combining 5 bias frames from each observation night with 
gbias, and normalized fiat frames were generated using gsf lat. All science images and our 
spectral standard (LTT9239) were then bias-subtracted, fiat-fielded, and overscan subtracted 
with gsreduce. Wavelength calibration was performed using the task gswavelength and 
the CuAr lamp spectra (also processed with gsreduce) taken at the time of observation. 
The wavelength calibration from the arc spectrum was applied to the 2D science spectra, 
which were then rectified using the task gstransf orm. The sensitivity function is determined 
by running the spectral standard through gsstandard, which is then applied to the galaxy 
spectra with gscalibrate. 

We do not use gsskysub, and instead apply gsextract immediately. This choice was 
made because gsextract uses the routine apall, which includes a sky subtraction proce- 



^IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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dure. If one runs gsskysub and then gsextract, the resulting variance spectrum incorrectly 
contains low uncertainty for regions that originally featured heavy sky contamination, as the 
sky-lines have been removed in the previous step. We also found that the Gemini tasks do 
not correctly propagate the VAR and DQ planes. In gsextract, the input DQ plane is not 
used, and the output DQ plane is simply set to be equal to (i.e., 'good') everywhere. The 
VAR plane output is the square of the uncertainty spectrum determined by the call to apall; 
the VAR plane that had been propagated up to this point is simply discarded. We still use 
this spectrum to describe the uncertainties in our flux measurements as it is the same as we 
would have had, if we had performed the reduction procedure outside of the Gemini suite of 
packages. 



3. Emission-Line Analysis 

Deriving the gas-phase oxygen abundance and the sSFR requires accurately measuring 
the flux in specific emission lines of the SNe la host galaxies; throughout this paper we 
will refer to the gas-phase oxygen abundance as 'metallicity', unless the context requires 
more clarity. Our procedure for obtaining emission-line spectra from observed data closely 
follows that in T04, in which a model of the continuum flux is subtracted from the observed 
spectrum. This corrects for absorption features superimposed on emission lines. We perform 
this analysis on all of our data, which totals 208 spectra from 138 host galaxies. 

We begiij^ by correct ing our spectra for extinction in the Milky Way along the line of 



sight (jSchlegel et al.lll998l ). and then mask the spectra in the region about the Balmer lines 
(Ha, H/3, H7, m) and the forbidden lines ([O II] AA3726, 3729, [O III] AA4959, 5007, [N II] 
AA6548, 6584, [S II] AA6717, 6731). We match this continuum spectrum with model spectral 
ene rgy distributions (SEDs) drawn from a grid of Composite Stellar Populations (CSPs) 
by Bruzual fc Chariot (j2003 . hereafter BC03), the parameters of which can be found in 



Table [H The model spectra are redshifted to the observed host-2;, and the fluxes are convolved 
with a Gaussian to mimic both the instrumental dispersion and the host's intrinsic velocity 
dispersion (also listed in Table [1]). The best-fit CSP model from our grid is determined for 
the masked host spectrum through a least-squares analysis (see Figure [1]). The continuum- 
subtracted spectrum is then smoothed by subtracting off the average flux in a 200 pixel sliding 
window, with emission lines again masked from this average. Fluxes are then measured by 
simulatenously fitting Gaussians to each of the aforementioned 12 lines, where the set of 



■^For spectra obtained with the SDSS 2.5ni telescope, the wavelengths are first converted from their default 
values ('vacuum') to 'air'. 
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all Balmer lines and the set of all Forbidden lines are each constrained to have a single 
common line width and velocity offset, resulting in 16 free parameters for 12 lines. Making 
this assumption strengthens the detection of weaker lines. 



60 



I I 

Host Galaxy of SN2006fw (z = 0. 08358) 




Ha 



5230 5260 5290 5230 5260 5290 5230 5260 5290 
X (Angstroms) 
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7000 
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Fig. 1. — : Spectrum of the host galaxy {z = 0.08358, nir = 18.4) of SN 2006fw, obtained 
by BOSS (Olmstead et al., in p rep). Spectrurn is sho wn in the observed frame, with Milky 
Way extinction correction from ISchlegel et al.l (119981 ). Key emission lines used throughout 
this paper are labelled. The signal-to-noise ratio (S/N) of this spectrum at ~ 5000A is 13; 
note, however, that the S/N of the emission lines are typically much higher than that of 
the continuum. INSET: An example of the continuum subtraction procedure outlined in 
Section [3J On the left is an enlargement of the full spectrum plotted below around the HP 
line; the middle figure shows the best-fit model for the continuum; and the right shows the 
continuum-subtracted spectrum, which we refer to as the 'emission- line' spectrum, and from 
which we measure fluxes. 



If the spectrum of a galaxy is dominated by an active galactic nucleus (AGN) instead of 
star formation, then emission-line ratios do not accurately reflect the gas-phase metallicity. 
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For this reason we use a BPT diagram ([O IIl]/H/3 vs. [N II]/Hq;; [Baldwin et al.lll98lf ) to 
identify and remove galaxies with AGN contamination from our host sample . We use the 



divisio n in this plane between star-forming galaxies a nd AGN calibra t ed by iKewley et al. 



fl200ll ). as well as the more conservative division by iKauffman et al. fl2003h: galaxies for 



which the two calibrations disagree are labelled 'composite' (similarly to iBrinchmann et al. 



20041 ). We perform our subsequent analysis both with and without the 'composite' hosts. 
We adopt the T04 S/N requirements; the Ha, H/3, and [N II] A6584 lines are required to be 
detected at > 5a, and [O III] > 3cr. To avoid the removal of high-metallicity hosts (which 
have a large log [N Il]/[0 III]), spectra with [O III] < 3a but log(Ha/[N II]) < -0.4 are 
retained in our sample. 



The remaining galaxy spectra are then corrected for host-ga^ 
that the intrinsic Ha/H/3 ratio follows Case-B Recombinatio n 



axy extinction by assuming 



Osterbrockl Il989l ) . and any 



deviations are due to extinction following the Cardelli et al. ( 19891 ) reddening curve, with 
Rv = 3.1. For the few cases where the observed Ha/H/3 ratio is bluer than the assumed 
intrinsic value, we make the assumption that Ay ~ 0. 

The spectra from Gemini are treated differently throughout this paper, as these only 
extend to ^ 66OOA in the observed frame. When the emission-line fluxes are fit we perform 
the same process as previously detailed, except we do not include constraints on Ha, [N II], or 
[S 11] (the Gemini hosts are all at 0.09 < z < 0.16). The lack of observed Ha and [N II] means 
we cannot use the BPT diagram to assess whether the spectrum is AGN or not. While visual 
inspection shows that none of the emission lines are broad, AGN contamination cannot be 
completely ruled out. For host-galaxy reddening corrections, the H7/H/3 ratio is used in 
place of Ha/H/3. We require H7 to have S/N > 10 for this correction, which holds true for 
most of our spectra. 

We summarize the statistics of our sample in Table [H We start with all 208 SN la 
host galaxy spectra, and then apply our cuts on the observed SN properties and on the 
host-galaxy spectral properties. The total number of spectra that pass all of our cuts is 
greater than number of SN host galaxies in our final sample due to multiple observations of 
the same host. 



3.1. Metallicity 



There exist many attempts at calibrating emission-line ratios to me tallicity, both em- 
pirica lly derived and based on photoionization models (for a summary, see iKewley Sz Ellison 
20081 . hereafter KE08). While the absolute calibrations of models rarely agree, relative dif- 
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Table 1. BC03 Composite Stellar Population Models 



Parameter Values 

IMF Chabricr 
[Z/H] 0.004, 0.008, 0.02, 0.05 

T [Gyr]^ 0, 0.1, 1.0, 4.0 

Dust'' 0, 0.25, 1.0, 2.0 

(Tint [km s-i] 50, 100, 150, 200, 250, 300, 
350, 400 

Age [Gyr] 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 

0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 
1.6, 1.8, 2.0, 2.5, 3.0, 3.5, 
4.0, 4.5, 5.0, 6.0, 7.0, 8.0, 
9.0, 10.0, 11.0, 12.0, 13.0, 14.0 



''Exponential Star Formation History, where 
SFR(t) ~ e-*/'^. '0' denotes a single burst (SSP). 

^ K-band optical depth for stars younger then 10 
Myr. 



Table 2. Breakdown of Spectral Host Observations 



Telescope 


Total 


Light Curve 


S/N 


AGN^' 




Spectra 


Cuts 


Cuts 


Cuts 


SDSS 


123 


59 


43 


28/10 


BOSS 


54 


31 


21 


17/2 


Gemini 


16 


12 


9 


9b 


NTT 


8 


6 


3 


3/0 


APO 


7 


7 


2 


2/0 



Note. — Light Curve cuts (Section [2} remove host- 
galaxy spectra from our sample for wrhich the SN doesn't 
have a distance modulus measurement from SALT2. Signal- 
to-noise and AGN cuts (Section [Sj remove spectra which 
cannot provide a metallicity or SFR measurement. Num- 
bers listed under the selection cuts are the quantity that 
passes each cut as well as cuts in the columns that precede 
it. 

^Spectra passing the AGN cuts are grouped as 
'SFG/Composite'; see Section |3] 

''As our observations with Gemini do not include Ha, 
we do not make any cuts for AGN activity. 
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ferences in metallicities are usually consistent - two different methods are likely to agree that 
one galaxy is more metal-rich than another. However, the calibration offset precludes mak- 
ing measurements in one system and directly comparing them to a measurement made via a 
different method. In KE08, fitting functions that transform metallicities computed in one sys- 
tem to those in another are determined for a varie ty of widely used emission-line techniques. 
We use the diagnostic of iKewley &: Dopital (120021 . hereafter K D02) throughout our analysis, 
except in the case of our Gemini observations, for which we use iKobulnicky fc Kewleyl ( 12004 
hereafter KK04). 

Galaxy metallicity is derived in KD02 through the ratio of [N II] A6584/[0 II], which 
is shown in their Figure 3 to be insensitive to variations in the ionization parameter at super- 
solar oxygen abundances. As discussed in Appendix A of KE08, for values of log [N Il]/[0 II] > 
— 1.2, the metallicity can be obtained as a root of the equation 

log([N Il]/[0 II]) = 1106.8660 - 532.15451Z + 96.373260^2 

- 7.8106123Z^ + 0.23928247Z^ (1) 

where Z = log([0/H]) + 12. The systematic uncertainty to these measurements is ~0.1 
dex (KE08). Fo r log [N Il]/[0 II] < -1.2, the average of the R23 calibrations by KK04 and 



McGaughl ( I1991I ) is used; this case applies to only 2 of our host galaxies. 



For the spectra taken at Gemini, we do not have observations of the [N II] lines, so 
we must use an i?23-based method (R23 = ([O II] AA3726, 3729 + [O IIl]AA4959, 5007)/H/3). 
We choose to use the KK04 calibration for metallicity derivations of these galaxies, which 
has a systematic uncertainty of ~0.15 dex (KE08). The R23 function is double-valued as a 
function of metallicity, with small values of R23 being representative of both low and high 
oxygen abundances. The KK04 method reaches a self-consistent solution for the ionization 
parameter and the metallicity by iteratively using the R23 line ratio and the parameter 
y = log ([O III] A5007/[O II]), the latter of which is necessary for constraining the ionization 
state. To break the degeneracy and determine to which of the two 'branches' of the double- 
valued function a galaxy belongs, the [N Il]/[0 II] ratio is typically used. Since [N II] is 
not observed in our Gemini spectra, we convert the measured -R23 value for each galaxy 
into metallicities from both branches, and explore in Section 14.11 the implications from all 
combinations of these metallicities. When we add the i?23-based metallicities from Gemini 
to the rest of our measurements, we use the KK04— 7-KD02 conversion function listed in 
Table 3 of KE08, which adds an RMS scatter of 0.05 dex (KE08). As described there, the 
transformation function is only valid for KK04-based metallicities with values 8.2 < Z < 9.2. 
All emission lines used in this analysis are labelled in Figure [U 



To truly measure the metallicity of the progenitor's environment, our spectra should 
be centered on the location of the SN. However, all of the spectra used in this paper were 
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not obtained in this manner; both the fiber spectra and sht spectra were centered on the 
galaxy core. The quantity that we are actually measurin g, then, is the avera ge metallicity 
from the integrated luminosity of a fraction of the host. iKewley et al.l ( 120051 ) showed that 
for an emission-line metallicity measurement to be representative of the global value, the 
spectrum should contain > 20% of the host-galaxy ^^-band fiux; KE08 showed that the g- 
band fiux fraction should be > 30% for galaxies with M > lO^^M©. Thus, the fraction 
of the total host-galaxy light which is in our spectra is an important quantity to consider. 
Spectra obtained as a part of the SDSS survey have fiber- and model-magnitudes associated 
with each galaxy, so we can easily compute the observed spectral fraction. As the BOSS 
fiber spectrum is only sampling the central 2" in diameter, we cannot use the SDSS fiber 
magnitude, which was obtained with a 3" fiber. The SDSS Catalog Archive Server (CAS) 
contains averaged surface brightness profiles within annuli of increasing radius. We use these 
brightness profiles to obtain the total apparent magnitude within the central 3 annuli, which 
span a diameter of ~ 2.05", and thus provides a good approximation to the BOSS fiber 
size. For the Gemini observations, inspection of the finder images clearly showed that all 
observed host galaxies had a majority of their light contained within the slit width. 



3.2. Star-Formation Rate 



The recent star-f ormation rate of a galaxy can be estimated by the fiux in the Ha nebular 
line (jKennicuttlll998l ) . whose prominence is due to absorption and re-emission of stellar light 
bluewards of the Lyman Limit. As only young, massive stars with M > IOMq contribute 
significant luminosity in this feature, the Ha diagnos tic is sensit i ve to star- formation in the 
past ~ lOVear s. We adjust the conv ersion factor in iKennicuttI (119981 ) for a Chabrier IMF 



(IChabrierll2003l ) by dividing it by 1.7 (jPozzettil 120071 ). which gives 



SFR [Mq yr' 



10 



-41.33 



L{Ha) [erg s 



(2) 



We conservatively assume a systematic uncertainty in log [SFR] of 0.2 based on lBrinchmann et al 



(120041 ). The if a- luminosity {L{Ha)) is derived from the measured fiux by assuming the best- 
fit cosmology detailed in Section [2l Equation [2] can only tell us the star-formation from the 
integrated light that comprises each spectrum; the SFR of the entire galaxy is also a function 
of the percentage of the total luminosity of the galaxy in each spectrum. For galaxies with 
a spectrum from SDSS or BOSS, we estimate the total SFR for the galaxy by scaling the 
quantity derived from Equation [2] by the inverse of the percentage of u-band light contained 



in the fiber (see iGilbank et al.ll2010l . Appendix A). This is obtained in a way analagous to 



that described for the (7-band light fraction in Section 13.11 
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We compute the sSFR of the host in two different ways. In the first method, we divide 
the SFR obtained using Equation [2] (the SFR within the fiber) by the mass derived from the 
least-squares fit of BC03 models to the continuum in Section [31 Each model is normalized 
to the luminosity from IMq of stars, so the multiplicative factor which minimizes the chi- 
squared of the fit for the data to each model (along with the distance modulus) provides a 
mass measurement. For every galaxy we make a cumulative distribution function (CDF) of 
the mass probability. We take the median mass value from this distribution as our estimate 
of the mass, with half the difference between the masses at 16% and 84% probability as our 
la error bar. As the resulting sSFR is obtained directly from our spectroscopic observations, 
it is denoted as sSFRspec. 

The second method we use to derive the sSFR is to take the SFR scaled to the entire host 
galaxy and divide by the full mass of the galaxy, using the mass as determined fr om fits of the 



observed UV, optical, and near-IR photometry to galaxy SED models, done in iGupta et al. 



(120111 ). The best-fit mass in their work is similarly defined as being the median of the CDF 
from their model fits. As this process requires the use of photometric measurements in 
obtaining both the mass and total SFR, we denote it sSFRphot- 



4. Results 

In Table [31 we list the derived oxygen abundances in 39 SN la host galaxies from useful 
spectra, excluding Gemini. All hosts in this table have a metallicity computed directly 
from the [N Il]/[0 II] ratio via the KD02 method. Specific star-formation rates from both 
methods detailed in the previous section are also listed. Although there are more emission 
lines required at a specified S/N for the derivation of the metallicity than for the SFR, if the 
host cannot be placed on a BPT diagram then we are unable to determine whether the lines 
come from a star-forming region or AGN. As such, we do not list any host spectra in Table [31 
with a SFR measurement and no metallicity. Where there are multiple spectra from SDSS, 
we give the weighted average for each quantity. We list measurements from different sources 
separately, as the derived quantities are dependent on the covering fraction of the spectrum. 
Also included in Table [31 are Hubble Residuals (HR) from SNANA using both mlcs2k2 and 
SALT2; the AGN contamination in the emission-line spectrum; and the fraction of the host's 
luminosity which is captured by the fiber, where applicable. 

A total of 16 host galaxies ai z < 0.15 were observed with Gemini (one of the galaxies 
had z = 0.1543). We are able to calculate an R23 value from the spectra in 11 of these 
galaxies, although two of these are host to SNe that fail the light-curve cuts for both salt2 
and mlcs2k2. For the two host spectra with the largest R23 values, the lower-branch metal- 
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licity derived from the fitting function of KK04 is higher than the upper-branch metalhcity, 
indicating that we are unable to constrain the metalhcity with this diagnostic. In Table H] we 
list the i?23 value; KK04 upper- and lower-branch metallicities; and the KD02-base metal- 
licities for both branches, using the conversion factors in KE08. These spectra are of strong 
emission-line galaxies, and as such our statistical uncertainties on R23 and all quantities 
derived thereof are small; the systematic uncertainties from the KK04 method and from the 
KE08 conversion function dominate our quoted errors. As in Table [Sj Hubble Residuals and 
sSFR measurements are listed in Table HI including for hosts which do not have metalhcity 
measurements. 

In th e analysis that fo llows we use the Markov Chain Monte Carlo (MCMC) package 



LINMIX flKelly et al.ll2007l ). which allows for uncertainties in both the x- and y-coordinates, 
to linearly fit for Hubble Residuals as a function of host metalhcity and sSFR. We use 
this package because our dependent variables are derived quantities, subject not only to 
measurement uncertainties but systematic uncertainties from the functions which relate the 
observable to the desired physical property. The linear regression coefficient s and their 



signifi cance are determined from the posterior distribution. Following the work of lKelly et al. 



( I2OIOI ). we cite both the percentage of runs from the posterior of the MCMC which have a 
slope with the opposite sign of the apparent correlation, and the significance at which the 
best- fit slope deviates from 0. 
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Table 3. Host Galaxy Metallicity and sSFR 



T A TT^ 




iOg[W/-rlJ-hlZ 


nrtMLCS 


n-ttSALT 


Sor -tXspcc 




fi'cictioii 




spectra,! 
Source 


2005OZ 


0.1298 


8.96(0.03) 


-0.29(0.11) 


-0.08(0.13) 


-11.98 


-11.52 


0.52 


1 


SDSS 


2005fr 


0.0902 


8.63(0.02) 


-0.07(0.09) 


-0.09(0.07) 


-10.42 


-9.77 


0.19 


1 


BOSS 


2005fv 


0.1182 


8.71(0.02) 


-0.01(0.09) 


-0.08(0.06) 


-10.44 


-9.91 


0.21 





SDSS 


2005fw 


0.1437 


8.60(0.03) 


0.18(0.07) 


0.04(0.06) 


-9.75 









NTT 


2005gp 


0.1266 


8.92(0.02) 


-0.20(0.12) 


-0.23(0.07) 


-10.56 


-9.75 


0.29 





SDSS 


2005hn 


0.1076 


8.85(0.02) 


0.09(0.11) 


0.03(0.07) 


-9.60 


-9.77 







APO 






8.91(0.01) 






-10.16 




0.47 





BOSS 


2005gb 


0.0866 


8.95(0.01) 


0.16(0.07) 


-0.04(0.06) 


-10.37 


-9.60 


0.17 


1 


SDSS 


2005ho 


0.0628 


8.62(0.01) 


0.04(0.08) 


-0.08(0.08) 


-9.81 


-9.26 


0.22 





SDSS 


2005hx 


0.1210 


f 


0.14(0.12) 


0.02(0.15) 


-9.96 


-9.05 


0.29 





BOSS 


2005if 


0.0671 


8.92(0.01) 


0.01(0.09) 


-0.12(0.09) 


-10.22 


-9.60 


0.16 





SDSS 


6213 


0.1094 


9.03(0.01) 




0.24(0.21) 


-10.39 




0.09 





SDSS 


2005ij 


0.1246 


8.83(0.01) 


-0.08(0.07) 


-0.12(0.05) 


-10.59 


-9.85 


0.20 





SDSS 


2005ir 


0.0764 


9.00(0.01) 


0.09(0.11) 


-0.04(0.07) 


-10.61 


-9.71 


0.14 





SDSS 


2005kp 


0.1178 


8.25(0.09) 


0.12(0.13) 


0.07(0.07) 


-9.50 


-8.85 


0.34 





BOSS 






8.37(0.03) 






-9.53 









NTT 


20051d 


0.1453 


8.72(0.01) 


0.08(0.12) 




-9.80 


-9.55 


0.22 





BOSS 


2006fc 


0.1217 


8.76(0.01) 




0.00(0.07) 


-10.21 


-9.53 


0.19 





BOSS 


2006fw 


0.0835 


8.81(0.01) 


0.12(0.08) 


0.01(0.06) 


-10.19 


-9.42 


0.37 





BOSS 


2006fy 


0.0827 


8.79(0.01) 


0.26(0.06) 


0.03(0.06) 


-10.10 


-9.63 


0.30 





SDSS 


2006fm 


0.1257 


8.85(0.01) 


0.12(0.08) 


0.06(0.06) 


-10.64 


-10.10 


0.10 





BOSS 


2006hl 


0.1482 


8.85(0.01) 


-0.04(0.08) 


-0.10(0.06) 


-10.60 


-9.98 


0.15 





BOSS 


2006hx 


0.0454 


9.10(0.02) 


-0.26(0.11) 


-0.22(0.10) 


-11.18 


-10.44 


0.22 


1 


SDSS 


2006kd 


0.1363 


9.06(0.02) 


0.29(0.08) 


0.28(0.06) 


-10.69 


-9.90 


0.29 





SDSS 






9.03(0.01) 






-10.69 




0.17 





BOSS 


15362 


0.1341 


8.64(0.03) 


-0.24(0.23) 


-0.02(0.27) 


-10.14 









APO 






8.86(0.01) 






-10.33 




0.57 





BOSS 


20061a 


0.1267 


8.44(0.02) 


0.49(0.09) 


0.45(0.06) 


-9.42 


-9.64 


0.53 





BOSS 


2006nc 


0.1240 


8.89(0.01) 


-0.03(0.13) 


-0.03(0.15) 


-10.79 


-9.74 


0.18 





BOSS 


2006nd 


0.1288 


8.95(0.01) 


-0.07(0.10) 


-0.07(0.06) 


-10.03 


-9.43 


0.28 





SDSS 


2007hx 


0.0798 


8.50(0.01) 


0.32(0.10) 


0.20(0.08) 


-10.43 


-10.01 


0.07 


1 


SDSS 






8.61(0.02) 






-10.11 




0.05 





BOSS 


2007jt 


0.1448 


8.96(0.01) 


0.04(0.08) 


-0.10(0.05) 


-10.18 


-9.82 


0.23 





SDSS 






9.01(0.01) 






-10.10 




0.14 





BOSS 


2007ju 


0.0636 


8.26(0.03) 


0.04(0.10) 


0.02(0.10) 


-10.33 


-9.73 


0.20 





BOSS 


2007jg 


0.0371 


8.79(0.01) 


0.35(0.13) 


0.16(0.12) 


-10.69 


-9.75 


0.10 





SDSS 


2007jd 


0.0727 


9.05(0.01) 


-0.01(0.13) 


-0.03(0.08) 


-10.68 


-9.94 


0.16 





SDSS 


20071g 


0.1104 


9.08(0.01) 


-0.17(0.08) 


-0.18(0.06) 


-11.12 


-10.75 


0.42 





BOSS 


20071o 


0.1384 


8.53(0.04) 


0.12(0.10) 


0.05(0.06) 


-10.26 


-9.76 


0.42 





BOSS 


20071c 


0.1150 


8.96(0.01) 


-0.04(0.09) 


-0.19(0.07) 


-10.75 


-10.03 


0.23 





SDSS 


2007ma 


0.1073 


8.94(0.01) 


0.03(0.08) 




-10.17 


-9.50 


0.29 





SDSS 


2007mh 


0.1278 


9.07(0.01) 


0.14(0.08) 


0.05(0.06) 


-10.70 


-10.02 


0.27 


1 


SDSS 


19048 


0.1368 


9.07(0.01) 


0.78(0.16) 


0.71(0.13) 


-10.26 




0.53 





SDSS 






9.08(0.01) 






-10.39 




0.43 


1 


BOSS 


2007mn 


0.0769 


8.83(0.01) 


0.20(0.07) 


0.14(0.06) 


-11.06 


-10.43 


0.16 


1 


SDSS 
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4.1. HR vs. Metallicity 

The (/-band fiber fractions for BOSS and SPSS sp ectra are listed in Table [3l which we 



use together with the guideline of iKewley et al.l ( 120051 ) to determine whether a given fiber 



spectrum is representative of a 'global' metalhcity. However, supposing that the metalhcity of 
the progenitor is correlated with the globally averaged value, a measurement from a region 
that is simply biased with respect to the average would still be useful. Including these 
measurements could add additional noise to our results, although this option would make 
our SFG sample almost 50% larger. In Figure |2] we examine the metallicity measurements 
from low fiber- fraction spectra. For SFG, these spectra on average measure slightly more 
metal-rich regions with a smaller dispersion than high fiber-fraction spectra, although this 
distinction disappears when 'Composites' are included. The low fiber-fraction sample does 
not include any extremely high metallicities when compared to the higher fiber-fraction 
sample. While we conclude that it is beneficial to use our full data set, for clarity we will 
quote our results from both a 'global' sample and 'full' sample in this section, where the 
former refers to ^f-band fraction > 0.20, and the latter sample includes all SDSS and BOSS 
measurements. For the slit-spectra in our data (NTT, APO, Gemini) the majority of the 
host was included in each observation, and these are included in both of our samples. 

For our observations using Gemini, we have upper- and lower-branch metallicity mea- 
surements, but lack the necessary data to discern the correct value. We resolve this issue 
by determining the best-fit linear relation between HRs and metallicity in a given sample 
for each possible combination of Gemini-observed host metallicities. The mean value of the 
significance over the whole set of possible metallicity combinations is taken, and we give 
limits based on the combinations which yield the lowest and highest significance. 

A correlation between metallicity and Hubble Residuals is measured with LINMIX at 
the 1.3 — 2.0(7 confidence level, depending on the combination of light-curve fitter and data 
sample used (see Table [5]). The trend is in the direction wher e lower metall i cities yield 



underlumino us SNe after light- curve corrections, as inferred in ISuUivan et al.l (120101 ) and 



predicted in iKasen et al.l ( 120091 ). We show the salt2 HR vs. metallicity plot for the 'full' 
sample of SFG in Figure |3l The hosts targeted by Gemini, being the faintest hosts within the 
redshift limit of our sample, unsurprisingly provide most of the weight in the low-met allicity 
region. The plot excludes the metallicity measurements from 2 of the 7 Gemini-observed 
host galaxies: those of SN 20061a and SN 20071o. Since both of these hosts also have BOSS 
spectra, which include direct observation of the Ha line, we used the derived quantities for 
these hosts in place of those taken from Gemini spectra. We note that the lower- and upper- 
branch Gemini metallicities of SN 20061a (8.36/8.39) agree with the BOSS value (8.44), 
and that the SN 20071o values, while different (8.16/8.79 for Gemini, 8.53 for BOSS), both 
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Table 3 — Continued 



lAU^' 


Redshift 


log[0/H]+12b 


HRmlcs" 


HRsalt" 


sSFRspcc^ 


sSFRphot'' 


(j-band 


AGN<= 


Spectral 


Name 














fraction 




Source 


2007OU 


0.1132 


8.91(0.01) 


0.48(0.10) 


0.45(0.08) 


-9.99 


-9.44 


0.23 





SDSS 


2007pd 


0.1399 


8.63(0.01) 


0.10(0.10) 


0.11(0.07) 


-9.95 


-9.97 


0.18 





BOSS 



''SNe in our Phot la sample do not have an lAU designation, and instead the internal SDSS-Supernova Survey candidate 
ID is given. 

''Metallicty values listed here arc in units of 12H-log [O/H]. For host galaxies where there is more than one spectrum from 
SDSS-I, we give a weighted average of our derived metallicities. Errors listed here are statistical uncertainties; a systematic 
uncertainty of 0.1 dex is added in quadrature for all analyses. 

'^An additiona l uncertainty of 0.14 magnitudes is added in quadrature to the listed uncertainty in all analyses (e.g., 
iKellv et allbOld) : this quantity is the intrinsic uncertainty in the SN la sample, required for the reduced of the best- 
fit cosmology to be 1. 

'*logsSFR[yr^^], derived from Ho flux. We adopt a systematic uncertainty of 0.2 dex for all SFR measurements. 

"'0' denote s a St ar-Forming Galaxy (SFG), '1' a Composite. Both are classified as star-forming by the criteria of 
iKewlev et al. 1 1 I2OOII) . but the latter are classified as AGN by that of iKauffman et al! 

^2005hx has log [N II] /[O II] < —1.2, but the R23 value is too high to constrain metallicity, as is also the case for two of our 
Gemini-based spectra. 



Table 4. Host Galaxy Metallicity and sSFR from Gemini Observations 



lAU Name 


Redshift^' 


log(fl23) 


•7 b 
^KK04 


lower branch 


ZkD02 
upper branch 


HRmlcs 


HRsALT 


sSFRspcc*^ 


2005hr 


0.1163 


0.908(0.006) 


8.45/8.51 


8.34(0.15) 


8.39(0.15) 


0.17(0.07) 


0.13(0.06) 


-9.22 


2005hx 


0.1210 


0.763(0.007) 


8.23/8.77 


8.14(0.13) 


8.66(0.17) 


0.14(0.13) 


0.02(0.15) 


-9.96 


2006jh 


0.1249 


0.737(0.003) 


8.13/8.81 




8.72(0.17) 


0.21(0.08) 


0.17(0.06) 


-9.82 


2006iz 


0.1363 


1.021(0.005) 








0.23(0.07) 


0.21(0.10) 


-9.24 


2006jq 


0.1276 


0.695(0.006) 


8.21/8.84 


8.13(0.13) 


8.74(0.17) 


0.11(0.07) 


0.04(0.05) 


-9.92 


2006kb 


0.1392 


0.810(0.009) 


8.35/8.68 


8.25(0.14) 


8.57(0.16) 


0.21(0.11) 


0.07(0.13) 


-9.68 


20061a 


0.1270 


0.903(0.005) 


8.47/8.51 


8.36(0.15) 


8.39(0.15) 


0.49(0.09) 


0.45(0.06) 


-9.42 


2OO7I0 


0.1386 


0.679(0.017) 


8.25/8.85 


8.16(0.14) 


8.76(0.18) 


0.12(0.10) 


0.05(0.06) 


-9.95 


20071t 


0.1140 


0.971(0.010) 








0.30(0.13) 


0.23(0.15) 


-8.86 



Note. — Metallicity values listed here are in units of 12-f log [O/H] 

''AH Hosts listed here have a statistical uncertainty in redshift of 0.0005 

''Metallicities derived assuming the lower/upper branches of the R23 relation. 

'^The lower limit of applicability of the KE08 transformation from KK04 to KD02 is 
12+log[O/H]KK04 = 8.2, which gives an upper limit to any lower value, in the KD02 system, of 8.11. 

''lo gsSFRfyr"-'-], derived from H/3 flux; we assume that the intrinsic Ha/H/3 ratio follows Case-B Recombination llOsterbrockl 
19891) . We adopt a systematic uncertainty of 0.2 dex for all SFR measurements. 
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Fig. 2. — : Metallicity measurements from SDSS-I and BOSS spectra as a function of host 
g'-band luminosity captured within the fiber. Blue circle s denote SFG, Red c ircles are Com- 
posites. The dashed line shows the division employed by iKewley et al.l (120051) to distinguish 
between spectra which represent a global average and those which do not. As expected, 
the average metallicity in low fiber-fraction spectra of SFG is higher than that in higher 
fiber-fraction spectra (black squares), though the difference is small (0.072 ± 0.039 dex). If 
the 'Composite' spectra are included, the average is the same to better than 0.01 dex. 
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suggest a lower metallicity than the median from our sample. 

The low-metallicity side of Figure E] contains almost exclusively hosts with positive HRs. 
We divide our samples into two metallicity groups about the line log[0/H] + 12 = 8.80, and 
compute the weighted average Hubble Residual and the error in the mean for each bin. We 
choose this value for o ur metallicity divisio n over a more natural choice of solar-metallicity 



(log[0/H] + 12 = 8.69; lAsplund et al.ll2009l ) for two reasons. First, our value better splits the 
two samples into roughly equal sized bins, yielding 16 (low-Z) and 18 (high-Z) SFG for the 
'full' sample, and 13 (low-Z) and 11 (high-Z) SFG for the 'global' sample. Second, we can 
only include Gemini-based measurements if both possible i?23-based metallicities fall within 
the same bin; a solar-metallicity division would leave out two host galaxies whose branches 
straddle the dividing line, while our choice retains all Gemini measurements. As can be 
seen in Figure El the HRs differ significantly between the two bins, with a weighted-average 
of —0.032 ± 0.016 in the high-metallicity bin and 0.082 ± 0.017 in the low-metalhcity bin. 
There is a 4.9a significance (0.114 ± 0.023) in the difference between the mean salt2 HR 
in SFG hosts in our 'full' sample between low- and high-metallicity environments. Similar 
significance is found if the smaller 'global' sample is used in place of the 'full' (0.116 ±0.027); 
if Composites are added to the SFG (0.091 ± 0.021); or if mlcs2k2 is used instead of salt2 
(0.132 ±0.031). 



In contrast with previous studies (IKelly et al.ll2010l : ISullivan et al.ll2010l : iLampeitl et al. 



20101 ). we do not include the intrinsic scatter in the Hubble Diagram (here ~ 0.14 mag) 
in our HR uncertainties. The intrinsic scatter is the quantity which must be added into 
the uncertainties of the measured distance moduli so that the reduced of the best-fit 
cosmology equals one. As such, it represents a physical process unaccounted for in our light- 
curve fitters that affects the overall luminosity of the SNe la. When attempting to explore 
the origin of this scatter, it is the measured difference between the distance modulus from 
each SN la, derived from known physical relations, and the distance modulus at its redshift 
for a given cosmology, that are of interest. Adding in an extra source of uncertainty in the 
measured Hubble Residual is equivalent to diluting the confidence in which the light-curve 
fitter is returning an over- or under-luminous SN. 



4.2. HR vs. sSFR 

We examined the correlation in actively SFG of both galaxy-averaged quantities, SFRphot 
and sSFRphot, with Hubble Residuals. We find no correlation between HR and SFRphot; the 
best-fit slope deviates from zero at < la for both the SFG-only sample and the SFG±Composite 
sample (as defined in Section [3]), whether we use mlcs2k2 or salt2 (See Table E]). This 
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Fig. 3. — : Hubble Residuals as a function of emission- line metallicity, measured using the 
KD02 method described in Section [31 Blue circles are SFG in our 'full' sample, excluding 
Gemini sources. The lower- and upper-branch metallicities from Gemini spectra are green 
and red squares, respectively. The plotted error bars include the systematic uncertainty 
of 0.1 dex added in quadrature to log[0/H] measurements (KE08). The orange diamonds 
represent the weighted averages of the upper- and lower-metallicity bins, with the division 
at 12 + log[0/H] = 8.80 set so that the bins contain approximately the same amount of 
data (18 and 16 hosts, respectively). We do not fold into our HR uncertainties the 'intrinsic 
scatter' of ~ 0.14 mag used to bring the reduced of the best-fit cosmology down to 1; see 
Section 14.11 for a discussion. 
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result is not unexpected; a given SFR could denote either a massive galaxy with a low sSFR 
or a small galaxy with a high sSFR. Thus it seems likely that a global SFR measurement, 
like the one used here, would sh ow no correlation w ith HR. Turning our attention to the 
sSFR, we use galaxy masses from I Gupta et al.l (120 111 ) to compute the HR-sSFRphot relation, 
shown in Figure HI The evidence for a correlation using LINMIX is slightly more significant 
here than for SFRphot but the relation is poorly described by our linear fit, as the derived 
sSFRs display both a tight distribution and high scatter. 
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Fig. 4. — : SALT2 HRs as a function of the photometrically defined specific star-formation 
rate, sSFRphot- Blue circles are emission-line hosts classified as star-forming galaxies (SFGs); 
red squares are Composites. Neither the best-fit line for the SFG-only sample (solid line) 
nor the SFG+Composite sample (dashed line) are good fits to the data, though they hint at 
a general trend of overluminous SNe being preferentially found in galaxies with a low sSFR. 



The results of our study of the correlations between the sSFRspoc and HR are also found 
in Table |6l As before, we show our results using both mlcs2k2 and salt2 HRs, and split 
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our samples into SFG-only and SFG+Composite. We show in Figure [5] the salt2 HRs as a 
function of our derived sSFRspec values. We find a correlation between HR and sSFRgpec at 
a confidence greater than 3cr irrespective of our choice of light-curve fitter, in the sense that 
overluminous SNe la after corrections tend to reside in hosts with lower sSFR. It is worth 
noting that only a fraction of the total intrinsic scatter in SN la luminosity is explained by 
this correlation. Subtracting the best-fit line in Figure [5] from the HRs of the SFGs, the 
intrinsic scatter is reduced from 0.14 to 0.11 mag. 




-12 -11 -10 -9 

Log{sSFR3p,, (Ha) [yr"^]} 

Fig. 5. — : salt2 HRs as a function of the spectroscopically defined specific star- format ion 
rate, sSFRgpec- Blue circles are hosts classified as star-forming galaxies (SFGs); red squares 
are Composites. The best-fit line for the SFG-only sample is shown. SNe la which are 
underluminous after light-curve corrections are preferentially found in regions of high sSFR. 
If divided into equally populated bins (green diamonds), the bin of lower sSFR has a typical 
HR consistent with that of passive hosts (brown arrow), the upper limit of which is show 
here. 

In the case of sSFRphot, the property that we are computing - a galaxy- averaged sSFR 
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- should be directly comparable from one host to another. However, the sSFRgpec mea- 
surements are based on observations that measure different proportions of the host's total 
luminosity, and as such are not obviously directly comparable. To explore what type of 
systematic effect this might create, we show in Figure E] our derived sSFRgpcc as a function 
of u-band fiber fraction. It is clear from this figure that our measured sSFRgpcc is not simply 
a function of the fixed physical fiber size from SDSS and BOSS. Thus the relationship we 
find between sSFRspec and HR cannot be dismissed as an artifact of the varying physical 
scales upon which our measurements are based. We do not correlate our measured SFRgpec 
with HR because the measured SFR, being a cumulative quantity, will depend strongly on 
the percentage of a galaxy's physical size enclosed in a given observation; this is why we 
previously chose to examine the global quantity SFRphot- 

Our results for sSFRgpec presented in Figure \5\ show a rapidly evolving relationship in 
sSFR, while sSFRphot does not. This result lends itself to two questions. First, why do the two 
differ? The sSFRphot measurement requires an extrapolation of the SFR, assuming that the 
SFR scales with the u-band magnitude. However, the lower the SFR in a galaxy, the higher 
the contribution to the flux in near-UV bands from older stars. Thus the approximation we 
make would tend to overpredict the SFR in galaxies with less star-formation, leading to the 
type of narrow distribution of sSFRs that appear in Figure |H 

The second question is whether the slope of the correlation between sSFRgpec and HRs 
is compatible with previ ous works. Figure O could be seen as being in contrast to the results 



of iLampeitl et al.l (120 lOf ) which, using photometric host-galaxy properties of SNe la from the 
SDSS-SNS, found a difference in HRs between SNe la in passive and SFG of ~ 0.1 magnitude. 
However, what we are measuring is the correlation wit h HR within the group of all SFG; 



a more accurate comparison with iLampeitl et al.l (120101 ) would be the mean HR of all SNe 



la in SFGs (plotted in Figure compared with that of SNe la found in passive galaxies. 
We define passive galaxies as those hosts where the Ha line is observed with S/N < 10, and 
which fail at least one of the emission-line cuts described in Section [31 For the 15 passive 
host galaxies with a 2a upper limit on sSFR below that of the lowest sSFR in Figure l5l we 
find an average HR of -0.037 ± 0.024. Compared to the average HR of 0.031 ib . 012 i n 



our SFG sample, this is a difference of ~ 2.5a and in agreement with iLampeitl et al.l ( 120101 ). 
Additionally, if we split the SFGs into equally populated upper and lower bins, we find an 
average HR of 0.111 ± 0.018 in the higher sSFR bin and -0.030 ± 0.017 in the lower sSFR 
one, the latter of which is consistent with what we find in passive galaxies. The similar HRs 
we find in both passive hosts and our lower bin of sSFR demonstrates that the linear fit 
plotted in Figure \5\ while a good representation of emission-line galaxies, does not capture 
the flattening of the trend at low sSFRs. 



Table 5. Significance of HR correlations with Gas-Phase Metallicity 



Sample SALT2 mlcs2k2 

'Global'iSFG 7.9%(1.3lgJo-) 5.8%(1.5lE!jo-) 

'Globar:SFG+Composites 8.4%(1.3tgj(r) 5.7%(1.5l|^j(r) 

'Fuir:SFG 2.0%(1.9lJ5;^o-) 2.8%(1.8l JJjo") 

'Full':SFG+Compositcs 2.0%{2.0l 0:30-) 4.0%(1.7to jcr) 



Note. — Entries are the percent of slopes in the posterior distri- 
bution of LINMIX with sign opposite that of the best-fit value, with 
the mean significance of the deviation from 0, based on all possible 
combinations of Gemini-based metallicities, in parentheses. The 
upper and lower limits on this significance reflect the strongest and 
weakest correlations for any single Gemini metallicity combination. 
Note that the significance of these slopes is lower in all cases from 
the significance at which HRs in the lower-metallicity bin differ 
from those in the higher-metallicity bin (see Section 14.111 . 



Table 6. Significance of HR correlations with SFR and sSFR 



Sample 


sSFRspcc 


sSFRphot 


SFR 


phot 




SALT2 MLCS2k2 


SALT2 mlcs2k2 


SALT2 


mlcs2k2 


SFG 


0.1%(3.1o-) < 0.1%(3.7(t) 


12.4%(1.1(t) 5.1%(1.6<t) 


19.9%(0.8<t) 


12.9%{l.lcr) 


SFG-I-Composites 


0.1%(3.2o-) < 0.1%(4.4cr) 


11.4%(1.2o-) 0.6%(2.6o-) 


21.2%(0.8o-) 


47.4%(<0.1o-) 



Note. — Entries are the percent of slopes in the posterior distribution of LINMIX with sign opposite that of the 
best-fit value, with the significance of the deviation from of the best fit in parentheses. The sample sizes differ 
between SFRpij^ti sSFRpij^t, and sSFRspoc- For example, 'SFG' for sSFRpiiot requires a photometric host mass 
and is limited to fiber-based spectra only, unlike sSFRspec- 
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Fiber Fraction of u — bond Host Luminosity 



Fig. 6. — : Specific star-formation rate measured within an SDSS-I or BOSS fiber, as de- 
scribed in Section 13.2^ as a function of tlie fraction of total n-band flux from the galaxy 
captured within the fiber. The absence of a correlation between sSFRgpec and fiber-fraction 
suggests that an apparent dependence of Hubble Residuals on sSFRgpcc is not a function of 
our measurements originating from different physical scales. 
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4.3. Completeness 

We have obtained spectra for nearly all of the host galaxies of SNe la at redshifts 
z < 0.15, covering the full range of host-galaxy absolute magnitudes. Combined with the 
completeness of the SDSS-SNS at these redshifts, we expect our study to be nearly free of 
selection bias with respect to both SN and host-galaxy parameters. However, having observed 
spectra for all of host galaxies does not mean that we are able to derive useful quantities 
from them all. A given galaxy spectrum can fail our emission-line S/N cuts for two reasons: 
either it is an elliptical galaxy and lacks these lines, or it is a low-luminosity galaxy. Only 
the latter group is important, as a 'complete' sample for our analysis by definition would 
exclude non-SFG. We lose some low-luminosity host-galaxies due to our data cuts, but the 
Gemini spectra provide us with a high S/N sample of intrinsically faint hosts to prevent our 
final analysed sample from being skewed. 

Figure [7] shows the cumulative distribution function (CDF) of SNe la host-galaxy r-band 
absolute magnitudes for all Spec la and Phot la SNe at 2; < 0.15. All gala xy magnitudes 



were obtained from the SDSS CAS, then de-reddened us ing the dust maps of ISchlegel et al. 



( 1l998l ) and k-corrected using the code kcorrect v4 . 1 . 4 (IBlanton fc Roweisll2007l ). Included 
in the plot is the host-galaxy CDF for all SNe la in this redshift range that pass the selection 
requirements for having a distance modulus measurement; it is clear that this cut introduces 
no bias into the host magnitude distribution. Also plotted is the sample which we draw our 
analysis from, the host galaxies that pass our S/N and AGN cuts (see Section [3]). We have 
far fewer bright galaxies, where < —21.5, but this is primarily due to the fact that many 
of these galaxies are passive. At the low-luminosity end, we are missing useful spectra for the 
faintest 5% of host-galaxies. This includes 4 of the 5 faintest hosts, 3 of which we observed 
with Gemini but were still of too low S/N for a metallicity measurement. A Kolmogorov- 
Smirnov test shows a 23% chance that the final host-galaxy magnitude distribution which we 
use in our analysis is drawn from the same sample as the sample which passes all SNANA 
cuts. Considering that we know of and expect the differences that do exist between the 
two distributions, it is reasonable to state that the host-galaxy properties in this paper 
approximate that of an unbiased sample. 



5. Discussion 



We have examined correlations between spectroscopic host-galax y properties in a nea rly 
unbiased sample with SN la Hubble Residuals for the first time fs eelKonishi et al.ll201ll for 
a recent analysis of a larger but biased sample with SDSS SNe; iGupta et al.l 1201 ll for an 



unbiased sample with multi- wavelength host-galaxy properties, also using SDSS SNe). We 
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Fig. 7. — : Cumulative distribution function in dereddened r-band absolute magnitude for 
host galaxies of 1) All SNe la, 2) SNe la that have distance modulus measurements /isN 
based on salt2, and 3) SNe la with both a /isn and a measurement of the host gas-phase 
oxygen abundance, based on emission-line diagnostics. All samples are for z < 0.15. The 
Kolmogorov-Smirnov test does not disfavor sample (3) - from which our final anaylsis is done 
- being drawn from the same distribution as sample (2). The differences which do appear 
between the samples are explained in Section 14.31 
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find evidence at > 4a of liiglier metallicity galaxies liosting overluminous SNe la (after light- 
curve corrections are applied) when compared to lower metallicity hosts, with a difference 
in the mean HR of ~ 0.1 magnitudes. However, the number of SNe la for which we have 
host-galaxy observations is insufficient to strongly constrain the nature of the relationship 
between the HR and metallicity, as evidenced by the much lower significance of a correlation 
returned by the LINMIX package. We also find at > 3a confidence evidence of galaxies with 
low sSFRs hosting overluminous SNe, with a difference of ~ 0.1 magnitudes between passive 
and star-forming hosts. 

The obs erved HR correlat ions are in the same direction as expected on theoretical 



grounds (e.g, iKasen et al.ll2009l . though larger in magnitude), so it is possible to take this 
as evidence of the overluminosity (after light-curve corrections) of SNe la from metal-rich 
progenitors. It is important to keep in mind, however, the difficulty in associating host- 
galaxy measurements with the properties of a SNe progenitors. The metallicity measurement 
used in this paper is indicative of the gas-phase oxygen abundance, which is most closely 
associated with episodes of recent star formation. While this should be strongly correlated 
with core-collapse SN progenitors, it is less clear how well this measurement relates to the 
older, lower mass progentiors of SNe la. In addition, whether we use the 'global' or 'full' 
sample in our analysis, the metallicity we are measuring is averaged over a large physical 
area, and in many cases is not derived from any galaxy light from the region where the SN 
occured. A useful future study would be centered not on the galaxy but the SN location 
itself; though given the long duration between the formation of the progenitor star and the 
detonation of the evolved white dwarf, it is uncertain that the location of the explosion is 
the region where the progenitor formed. 

Supposing that the SFR we measure is representative of the environment of the SN pro- 
genitor, what does a correlation between this quantity and SNe la HRs imply? Several studies 
have shown evidence of there being two populations of SNe la: 'prompt' a nd 'delayed', the 



monikers referring to the time elapsed sin c e the formation of th e progenitor flMannucci et al. 



20051 . 12OO6I : IScannapieco fc Bildstenll2005l : ISuUivan et al.ll2006l ). It is possible, then, that the 



galaxies we observe to have high sSFR are, on average, hosts of the 'prompt' population, 
while galaxies with a low sSFR are hosts to the 'delayed' population. However, the diagnos- 
tic which we use to measure SFR tells us about only the past ~ 10^ years, which is below 
the timescale on which even the 'prompt' SNe la progenitors transition from formation to 
explosion. 

It is also possible that the sSFR is simply acting as a tracer of metallicity, meaning that 
the correlation which we measure between HRs and sSFR is really one of HRs and progenitor 
metallicity. It is well known that for a given galaxy mass, gas-phase metallicity decreases 
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as SFR increases (see, eg. iMannucci et al.ll2010[ ). Thus, a positive HR-sSFR correlation is 



analogous to a negative HR-Z relationship, which are shown in Sections 14.11 and 14.21 Given 
the small number of host galaxies we analyze in this work, we are unable to determine 
whether our results are a result of this degeneracy, or whether there are two distinct trends 
of age and metallicity with HR. 

We note that the importance of the question, Ts the measured galaxy metallicity strongly 
related to the pro g enitor metallicity,' has led to it recently being addressed theoretically by 



Bravo fc BadenesI (120111 ). They find that the host metallicity does provide a good estimate 
of the SN metallicity in the case of an actively star-forming galaxy, which are the type of 
hosts that have been analyzed in this work. 

Given that we are using integrated galaxy spectra in an attempt to constrain properties 
of a single star - one which has since ceased to exist - it is inevitable that there be caveats 
and uncertainties regarding the precise relationship between the observable and the object 
of interest. But the effort is worthwhile, as the origin of the intrinsic error in the Hubble 
Diagram and any bias it may induce as a function of redshifts will contribute a significant 
systematic uncertainty that upcoming wide, deep surveys (DES, LSST) will face in the 
pursuit of constraining cosmology with SNe la. 
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